Universal local pair correlations of Lieb-Liniger bosons at quantum criticality 
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The one-dimensional Lieb-Liniger Bose gas is a prototypical many-body system featuring universal 
Tomonaga-Luttinger liquid (TLL) physics and free fermion quantum criticality. We analytically 
calculate finite temperature local pair correlations for this model at quantum criticality using the 
polylog function in the framework of the Yang- Yang thermodynamic equations. We show that the 
local pair correlation g^(0) has the universal value g' 2 '(0) « 2p/(ne) at quantum criticality, where 
n is the linear density and p is the pressure per unit length rescaled by the interaction energy 
e = rj^c 2 with interaction strength c. This suggests the possibility to test finite temperature local 
pair correlations for the TLL over the relativistic dispersion regime and to probe quantum criticality 
with the local correlations beyond the TLL phase. Furthermore, thermodynamic properties at 
high temperatures are obtained by both high temperature and virial expansion of the Yang- Yang 
thermodynamic equation. The virial expansion is more suitable for the weak interaction regime. 

PACS numbers: 03.75.Ss, 03.75.Hh, 02.30.1k, 05.30.Rt 



I. INTRODUCTION 

Advances in the trapping and cooling of atoms in op- 
tical waveguides have opened up exciting possibilities for 
testing theory in low-dimensional quantum systems. Ob- 
served results to date are seen to be in excellent agree- 
ment with results obtained using the mathematical meth- 
ods and analysis of exactly solved models [H-0I- These 
include the remarkable experimental realization of the 
Tonks-Girardeau gas [H-Q, the super Tonks- Girardeau 
gas 0], Yang- Yang thermodynamics on an atom chip 
8|, and the phase diagram of the attractive Fermi gas 
9(. Such exquisite tunability with tightly confined ul- 
tracold atoms provides unprecedented opportunities for 
improving our understanding of novel quantum phenom- 
ena such as quantum criticality, universal scaling the- 
ory, spin-charge separation and Tomonaga-Luttinger liq- 
uid (TLL) physics. 

The one-dimensionsal (ID) delta- function interacting 
Lieb-Liniger Bose gas [l| is a many-body system solved 
exactly by the Bethe ansatz hypothesis. It has had a 
tremendous impact as an archetypical system in quantum 
statistical mechanics [lOj]. It's i?-matrix provides 
the simplest solution of the Yang-Baxter equation [T^ |. 
Yang and Yang [l| showed that the thermodynamics 
of Lieb-Liniger bosons can be determined from the min- 
imisation of the Gibbs free energy subject to the Bethe 
ansatz equations [3, [l5j]. This thermodynamic Bethe 
ansatz method has been extended to a wide range of ID 
quantum many-body systems [lH ]. In particular, Yang- 
Yang thermodynamics is fundamental to the Y-system 
which has emerged as a ubiquitous integrable structure 
in mathematical physics [l6| . In the present context, it 



provides the framework to study thermodynamics, quan- 
tum criticality and TLL physics. The equation of state 
for Lieb-Linger bosons has been obtained [l7| analyti- 
cally for strong coupling and low temperature in terms 
of the polylog function. The expression for the equation 
of state in this regime enables the exploration of TLL 
physics and quantum criticality in this system. 

In a recent experiment [Isj . thermal fluctuations were 
studied in a highly elongated weakly interacting Bose gas 
at high temperatures, where the quantum fluctuations 
are strongly suppressed. In this regime, the measured 
thermal fluctuations are in good agreement with the ex- 
act Yang- Yang thermodynamics. However, at quantum 
criticality, quantum fluctuations are strongly enhanced. 
It is of particular interest to understand quantum cor- 
relations and fluctuations at quantum criticality. A fi- 
nite temperature quantum phase transition does not ex- 
ist in the ID Lieb-Liniger Bose gas. However, there is a 
critical point in the grand canonical ensemble, when the 
chemical potential /i c = 0, which separates the vacuum 
from a filled "Fermi sea" of particles at zero temperature. 
At finite temperatures, a TLL with relativistic disper- 
sions can be sustained in a region of the T — fj, plane. 
This implies that for temperatures below a crossover 
value T*, the low- lying excitations have a linear rela- 
tivistic dispersion relation. If the temperature exceeds 
this crossover value, the excitations involve free quasi- 
particles with non-relativistic dispersion. This crossover 
temperature can be identified from the breakdown of 
linear temperature-dependent entropy, see Fig. 1. In 
this phase diagram, quantum criticality is in the regime 
where t = ksT '/ '(^c 2 ) is small, but ksT > \fi — fi c \, or 
for the temperature below the degenerate temperature 
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FIG. 1: Quantum phase diagram of Lieb-Liniger bosons. The 
plot shows the entropy in the T — p plane obtained from 
the Yang- Yang thermodynamics. A crossover temperature 
separates the TLL critical phase with dynamic exponent z = 1 
and correlation length exponent v = 1 from the quantum 
critical regime (QR) with exponents z = 2 and v = 1/2. The 
other crossover temperature separates the quasi-classical gas 
from the quantum critical regime. 
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Here kg is the Boltzmann constant. 



The local two-body (pair) correlations have been ana- 
lytically calculated for a few limiting cases at T = [l9| 
and at finite temperatures from the Yang- Yang thermo- 
dynamic equations [20M2^ | . However, the local pair cor- 
relations have not been thoroughly investigated at quan- 
tum criticality. In principle, the Yang- Yang thermody- 
namic equations enable the numerical calculation of the 
local pair correlations in the whole parameter space. The 
local pair correlations are related to the universal con- 
tact which measures the probability of two bosons (or 
two fermions with opposite spins) at the same position 
25, 261. 



In the present paper, we calculate the local pair cor- 
relations of Lieb-Liniger bosons in analytic fashion from 
the Yang- Yang thermodynamic equations using the poly- 
log function in order to understand universal features 
of quantum fluctuations and critical phenomena in an 
archetypical quantum system. In addition, a high tem- 
perature expansion of the Yang- Yang thermodynamic 
equations is obtained in the strong coupling regime. This 
allows one to extract the contributions from quantum 
and thermal fluctuations to the classical Boltzmann gas 
at high temperatures. 



II. MODEL AND EQUATION OF STATE 

The Hamiltonian 
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of the Lieb-Liniger Bose gas [l[ describes N spinless 
bosons with mass m constrained by periodic boundary 
conditions on a line of length L. The contact interactions 
are governed by the effective ID coupling constant g\u = 
— 2h 2 /(mai]j) where a%D = (— a]_/2a s ) [1 — C (a s /a±)] 
is the ID scattering length with a± = y/2h/mu± and 
the numerical constant C f=s 1.4603 27]. The dimen- 
sionless interaction strength is defined by 7 = c/n with 
c = —2/aiu, where n — N/L is the linear density. The 
interaction strength can be controlled by tuning either 
or a s in experiments. 
In describing the thermodynamics of the model the key 
quantity is the dressed energy 



e(k)=Tln(p h (k)/p(k)) 
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which plays the role of excitation energy measured from 
the energy level e(fcp) = 0, where kp is the Fermi-like 
momentum. The thermodynamics of the model in equi- 
librium follows from the Yang- Yang equation [l3T ] 

/oo 
dqa 2 (k-q)ln(l+e-< q ^ T ) (3) 
-co 

where e°(k) — j— k 2 is the bare dispersion, /1 is the chem- 
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The pressure p(T) and the free energy F(T) are given 
in terms of the dressed energy by 
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F(T) = pn- 
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dk ln(l +e~ e W/ T ). (6) 



The pressure of the strongly coupled gas at finite tem- 
peratures has been obtained from ([3]) in the form [17j 
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where Li s is the standard polylog function and 
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Furthermore, the Yang- Yang equation ^ can be ex- 
panded in powers of the dimensionless temperature t = 
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k^T je with e = (^-c 2 ). Thus from Q the dimensionless 
pressure p = p/(ec) at finite temperatures follows as 
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with 



A = f L -±Lu 3 _(-e A °/ t ) + -^Lis(-e A V*) (i ) 



and 



A = A - -^Li^-e^*). 



(11) 



The result is essentially a high precision equation of 
state for Lieb-Liniger bosons at quantum criticality. We 
will verify that it is also valid for the high temperature 
regime as long as /cbT *C s. Recalling the phase diagram 
Fig. Q] the density and the compressibility can be cast 
into universal scaling forms [13, SSI 



n(T,fi) - n (T, fi) « U +1 
n(T,n)-Ko(T,ii) « ^ +1 ~^Q 
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near the quantum critical point fi c — 0. Here the dy- 
namic exponent z = 2 and the correlation length expo- 
nent v = 1/2 with the scaling functions given by 



T(x) = Lii(-e*) 
2v 7r 



Q(.x) 



2eV^F 



LU(-e") 



(14) 
(15) 



for f > |/z — /i c | in dimensionless units. The background 
density and compressibility in the vacuum are zero, i.e., 
no(t,n) = Ko(t,n) = 0. These analytical results provide 
insight into quantum fluctuations near the quantum crit- 
ical point. We will further demonstrate that the finite 
temperature local pair correlations shed light on quan- 
tum critical behaviour. 



III. LOCAL PAIR CORRELATIONS 

The general calculation of correlation functions in 
quantum many-body systems is a notoriously difficult 
problem, even more so at finite temperature [14] . Fortu- 
nately the two- particle local pair correlation g^(0) can 
be calculated relatively easily from the free energy. In 
the grand canonical description, the Hellmann-Feynman 
theorem gives 
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FIG. 2: Local pair correlation g^ 2 \0) at low temperatures 
in the t — fi plane. The right black-dashed line denotes the 
crossover temperature separating the TLL phase from the 
quantum critical regime (CR). The left black-dashed line sep- 
arates the vacuum (V) from the quantum critical regime. 
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FIG. 3: Local pair correlation <?' 2 '(0) vs chemical potential 
at quantum criticality. The correlation function </ 2 '(0) at 
different values of temperatures show a good agreement be- 
tween the analytical relation (|20jl and numerical result ob- 
tained from the TBA equation ((3j . 



where / = fj, — p/n is the free energy per particle. For 
constant n the local pair correlations require the calcu- 
lation of the derivatives dp,/dc and dp/dc, which can be 
done by iteration. Using the result ([7]) we find 



dp 
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where dA/dc can be calculated from (TTUl) by iteration. 
The calculation of d/i/dc is cumbersome. To this end, 
we first find the explicit form of the density n = dp/d/j, 
and then take the derivative with respective to c. After 
a lengthy iteration, we find 
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Finally, from the expressions (fTTj) and <jTHJ) , we obtain 
the local pair correlations in the form 
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where A is given by (|10| . 

We note that the local pair correlations satisfy a simple 
relation 
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where t is small. It turns out that this simple looking 
result not only holds for the quantum critical regime, but 
also for the TLL phase and the ground state at T = 0. In 
this sense it is a universal relation. The low temperature 
local correlation for the TLL phase follows from (|20p as 
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which coincides with the result given in 20-23]. 

Fig. [2] shows a plot of the local pair correlations ob- 
tained from (|19p in comparison with the numerical re- 
sult. The crossover temperatures, which separate the 
quantum critical regime from the TLL, are determined 
by the breakdown of the TLL with the local correlation 
pTjl . The crossover line is consistent with the phase dia- 
gram Fig. [TJ The TLL phase persists below the crossover 
temperatures, where both results (|2U)) and (|21l) coincide, 
see Fig. [2j This represents a smooth crossover from the 
rclativistic TLL regime to the regime governed by a non- 
rclativistic dispersion relation. The local correlations in- 
crease as the chemical potential becomes large and posi- 
tive because of the decrease of the interaction. The tem- 
perature enhances the local pair correlations, see Fig. [3] 
We see clearly that the local correlations (|2"0)l cover the 
quantum critical region, the TLL phase and ground state 
at T = 0. 

For nonzero temperature, the vacuum can be taken as a 
semi-classical gas regime, where the particle density n ~ 
_i e -H/T with thermal wavelength A" 1 = ^/mk B T/2Trh 2 , 
which is much smaller than the mean distance between 



two particles. The local pair correlation tends to zero. 
We see that the relation (|20l) between the local correla- 
tions and pressure holds in the physical regime as long as 
the temperature is below the degenerate temperature and 
T <C ^-c 2 . Near the critical point fx c = 0, the density 
is very low, therefore the interaction is strong and sits in 
the Tonks-Giraradau regime. At high temperatures, this 
crossover disappears due to suppression of the quantum 
fluctuations. At high temperatures, the chemical poten- 
tial become more negative. We consider this limit in the 
next section. 



IV. HIGH TEMPERATURE EXPANSION 

The recent measurements [18] on thermal fluctuations 
in the ID Bose gas were carried out in the weak coupling 
and high temperature regimes. The variance of atom 
number (SN 2 ) in a volume A can be evaluated as [l8| 



(<W 2 ) = Ak B T 
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where the local density n(fi, T) can be determined from 
the equation of state in the thermodynamic limit. Here 
we first consider the equation of state of this system at 
high temperatures and in strong coupling regimes. At 
high temperatures, the TBA equation (J3J can be ex- 
panded in an appropriate form. For convenience, we de- 
fine the function rj(k) — e e ( k ^ T and the inverse tempera- 
ture parameter a = h 2 c 2 /(2mT). The TBA equation Q 
can then be written 

^{k) = x k e^ dqa2{k - q)i < 1+ ^ 1(q) ) , (23) 

where Xk — e^^ e °^ k '^' T <C 1 at high temperatures. We 
carry out an expansion of this equation in powers of Z = 
e n/T <j; i as T -> oo. Collecting the first few terms, we 
have 

^(k) a x k (1 + BtZ + B 2 Z 2 + B 3 Z 3 ) + 0{Z b ) (24) 
where the coefficients Bi, for finite a, are given by 
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Using the expansion (|24[) gives the pressure ([3]) in the asymptotic form 



2nh 2 



T-^Z [1 + hZ + f 2 Z 2 + hZ 3 \ + 0(Z 5 ), 



(25) 



(26) 



with 
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where Erf(x) is the standard error function. 
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The TBA results reduce to several limiting cases. The 
pressure p c of the classical Boltzmann gas 
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follows in the limit T — > 00. The pressure of the ideal 
Fermi gas 
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is obtained in the limit c — > 00. Similarly the pressure of 
the ideal Bose gas 
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T 3 / 2 Li3(e^ T ) 



(30) 



follows in the limit c — > 0. 

In the high temperature limit, the chemical potential 
tends to negative infinity. In this region, far away from 
criticality, the pressure (0 given in terms of the polylog 
function and the high temperature expansion result (|2"6")l 
are highly accurate for a > 5, see Fig. |4j The polylog 
function result (0 gives a better fit with the numerical 
result obtained from the TBA © even for large values 
of fegT/fi. At high temperatures, although the thermal 



fluctuations dominate, the quantum statistics are still mi- 
croscopically significant. The quantum statistical effect 
is evidenced from changing the value of a, which can 
be controlled in current experiments [l8|. In the Tonks- 
Girardeau limit a — ¥ 00 the gas approaches free fermions. 
For small values of a, i.e., in the week coupling region, 
we see the pressure (|26[) obtained from the high temper- 
ature expansion is consistent with the numerical result 
at high temperatures, see Fig. But as expected, the 
result (|26|) is no longer accurate for a -C 1. 

In the weak coupling limit and at high temperatures it 
is more practical to consider a virial expansion with the 
TBA equation (0, i.e., 
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with Z = exp(/i/T). The first term in the rhs of this 
equation comes from free bosons and 
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contributes to the second virial coefficient. From the 
comparison shown in Fig. [5j we see clearly that the pres- 
sure obtained from the virial expansion gives a high pre- 
cision equation of state for a < 1. This result is valid for 
the experimental setting with a < 1 18]. 



6 



; 0.96 - 



I 0.94 
0.92- 
0. 




Numerical cc=5 

o High-T expansion a=5 

x Polylog a=5 
Numerical cc=15 

• High-T expansion 
Polylog a=15 

Classical gas 

Ideal fermi gas 



0=16 



0.2 



0.3 



0.4 



0.5 0.6 
-K T/u 



0.7 



FIG. 4: The ratio of the pressure p/p c vs ftsT/fi at high 
temperatures for a — 5 and 15. Here p c is the pressure of 
the classical Boltzmann gas. High accuracy of the pressure 
is seen for — oc < fcsT//x < 1. At high temperatures both 
pressure © obtained from the polylog formalism and the one 
(|26|) obtained from high temperature expansion are accurate 
by comparing with numerical result from the TBA 0. But 
the polylog formalism is also valid at low temperatures as long 
as a > 1. 
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FIG. 5: Comparison between different values for the pressure 
p/p c vs ksT '/ fi at high temperatures for a — 0.1, 0.2, 1.0. For 
small values of a, the pressure ()31[) obtained from the virial 
expansion is in good agreement with the result obtained from 
the numerical solution of the TBA equation ([3]). 



lations for the TLL over the whole relativistic dispersion 
regime and to test quantum criticality with the local cor- 
relations beyond the TLL phase. In the quantum critical 
regime, the thermal fluctuations strongly couple to the 
quantum fluctuations with universal free fermion z = 2 
quantum criticality. The local pair correlations provide 
insight into the microscopic quantum statistical effects 
at quantum criticality. We also derived the thermody- 
namic equations ((26)) and (f31"j) of the Lieb-Liniger gas at 
high temperatures using a high temperature expansion 
of the Yang- Yang thermodynamic equations. The effect 
of quantum statistics is microscopically significant even 
in the thermal fluctuation dominated high temperature 
regime. In particular, the equation of state given by (|3l"j) 
is high accurate for the weak coupling and high temper- 
ature regimes. 

Our analytical prediction (l2Tfl) for local pair correla- 
tions can be tested using current experimental techniques 
for preparing and detecting ID gases. An ensemble of 
parallel ID Bose gases can be prepared 0, 0] by loading 
ultracold Bose atoms into a two-dimensional (2D) optical 
lattice. The 2D optical lattices, an ensemble of parallel 
ID tubes, can be formed by superimposing two standing- 
wave lasers on the crossed dipole trap. The depth of the 
2D lattice must be sufficiently large to make the quan- 
tum tunnelling among these ID tubes negligible. That is, 
atoms in these ID tubes are almost all in the lowest trans- 
verse vibrational state. The photoassociation techniques 
for measuring local pair correlations in zero-temperature 
ID Bose gases could also be used to explore local pair 
correlations in non-zero temperature systems, such as the 
Lieb-Liniger Bose gas at quantum criticality studied in 
this paper. In particular, the universal relation between 
the local pair correlation and the pressure at quantum 
criticality could be explored following the experimental 
scheme for measuring the homogeneous contact of a uni- 
tary Fermi gas [32| . 



V. CONCLUSION AND DISCUSSION 

In conclusion, we have studied finite temperature local 
pair correlations of the Lieb-Liniger Bose gas at quan- 
tum criticality. A simple relation (|20|) between the local 
correlations and the pressure has been obtained in the 
framework of the Yang- Yang thermodynamic equations. 
This relation holds for both the quantum critical regime 
and the TLL phase at quantum criticality. It provides a 
simple way to probe finite temperature local pair corre- 
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